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. . Abstract 

We assume that dark matter is composed of scalar particles that form a Bose-Einstein condensate 

■ (BEC) at some point during the cosmic evolution. Afterwards, cold dark matter is in the form of a 

■ condensate and behaves slightly different from the standard dark matter component. We study the 
large scale perturbative dynamics of the BEC dark matter in a model where this component coexists 
with baryonic matter and cosmological constant. The perturbative dynamics is studied using neo- 
Newtonian cosmology (where the pressure is dynamically relevant for the homogeneous and isotropic 
background) which is assumed to be correct for small values of the sound speed. We show that BEC 

■ dark matter effects can be seen in the matter power spectrum if the mass of the condensate particle 
lies in the range 15meV < < 700meV leading to a small, but perceptible, excess of power at large 
scales. 

o 

^' 

Standard cosmology relies on the assumption that dark matter (DM) represents around 25% of the 
Q ■ cosmic energy budget and behaves as a pressureless component. This behavior allows the gravitational 

clustering of the DM particles (once they are decoupled from relativistic species at the early Universe) in 
^ ■ order to form dark halos that host galaxies. Such large scale structure supports the hypothesis that DM is 

composed of weakly interacting massive particles (WIMPs) though the real nature of this component 
is still unknown. However, the WIMP scenario faces several challenges at galactic scales [2]. 

The real nature of the dark matter is unknown, i.e. we do not know what it consists of. One possibility 
is the existence of dark scalar particles (i.e. a spin-0 bosons). In this case, the scalar field (p associated 
\ with such particles has its dynamics governed by the potential V{(t)) which encloses all the interactions 

■ of such system. Hence, if dark matter particles are bosons they can form, at some critical redshift Zcr, a 
O^l I Bose-Einstein condensate (BEC), i.e. ground state interacting bosons trapped by an external potential 

[3]. Such BEC stage of the dark matter was proposed in Refs. [H [5] and it has been widely studied in 
the literature ^ (see also [7]). 

In condensed matter physics the mean-field approximation is widely used. In such case, BEC systems 
are studied through the time-dependent generalized Gross-Pitaevskii equation [5] 



1 Introduction 



ih—^-^V^^ + V{r,t)^ + gm^)^, (1) 

where ^ = 5'(r, i) is the macroscopic wavefunction of the condensate, is the mass of the particle, 
V{r^ t) is the trapping potential. The non-linearity term reads 

g(|*|2) =wo|*|' + e|*|', (2) 

where the quadratic term accounts for the two-body interparticle interaction with ug = ATrh^la/m^. The 
term e is strength of the three-body interparticle interaction which is neglected in numerous works on 
BEC. The quartic term becomes important only in the case of higher densities [9]. In the standard 
approach one considers the case where e = 0. The free parameters are the boson-boson scattering length 
{la) and the mass of the particle (m^). 

The application of BEC physics in astrophysics has become much more clear in Ref. [10] where the 
authors have assumed an arbitrary non-linearity term with being described by the so-called Madelung 
representation [11], 

* = VpMJe^^^'-'*) (3) 

where p = is the density of the condensate. Hence, it comes out that BEC dark matter can be 
described in terms of a non-relativistic Newtonian fiuid with barotropic equation of state that obeys the 
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energy and momentum balances. The BEC pressure is given by |10j 

Pbe =g{p)p- j 9{p)dp. (4) 

In the standard approach of BEC systems g — uq I^*]^, the effective pressure of the condensate becomes 
Pbe P^- Within this approach remarkable results concerning galaxy curve rotation |10j . dynamics of 
galaxy clusters [12], the core-cusp problem [T3] and galactic vortices [3] were obtained (see also [T5]). 

An important issue that should be addressed is whether the background evolution and the large 
scale structure of the Universe can be affected by the existence of such BEC phase of the dark matter. 
In reference [TB] it was assumed that the BEC takes place via a first order phase transition once the 
temperature of the bose gas reaches the critical value Tcr- The cosmological parameters at the moment 
in which such event occurs (i.e. density of the bose gas, the temperature and redshift) were also established 
in terms of the free parameters of the BEC, namely the mass of the dark matter particle and la- 

The growth of BEC dark matter perturbations was also considered in references [171 UHl US] ■ The 
general conclusion is that assuming a positive scattering length (that means a repulsive self-interaction) 
the growth of BEC dark matter inhomogeneities is faster than the standard CDM case {la = 0). This 
would lead to an excess of power at small scales though the bottom-up structure formation scenario still 
works. In principle, such excess drives the large scale in a different way resulting in a different power 
spectrum P{k). As an example, ultra light bosonic DM models (m^ ^ lO^^^eV) are able to leave very 
small imprints into the acoustic peaks of the CMB [50]. However, it has been claimed that although a 
different small scale clustering the collapsed BEC structures (bosons stars, black holes, etc) behave as 
cold dark matter at large scales without modifying the large scale power spectrum [19| . 

The purpose of this paper is to investigate more carefully the consequences of a BEC dark matter 
cosmology for the matter power spectrum. Does a Bose-Einstein condensate leave imprints into the P(k)l 
Using the 2dFGRS power spectrum data [21] we will be able to investigate scales O.OlMpc""^ < kh~^ < 
0.185Mpc~^ that still belong to the linear perturbation theory. We will assume a cosmological background 
evolution as proposed in [16] . where BEC dark matter coexists with baryons and cosmological constant. 
We reproduce such model in section 2.1. In section 2.2 we introduce the cosmological perturbation theory 
using the neo-Newtonian equations. In order to compute the matter power spectrum, perturbations in 
both dark matter and baryons have to be taken into account. Indeed, the observable power spectrum 
is the statistical distribution of the visible matter (baryons) which has been driven by the dark matter 
gravitational field. This step has not yet been done for the BEC scenario. In section 3 we compare the 
BEC matter power spectrum with the data assuming different model parameters. We conclude in section 
4. Throughout this paper, we shall assume c = 1. 



2 Cosmological dynamics of Bose-Einstein condensates 
2.1 The background evolution 

After the BEC forms its effective pressure can assume a polytropic equation of state such as Pbe ~ p2e 
if one assumes an arbitrary non-linearity term . The exact value of 7 is defined by the non-linear 
contribution of the Gross-Pitaevskii equation which in its standard form leads to [lOj 

Pbe = —Pbe- (5 

ml 

The scattering length (la) and the mass {m^) of the dark matter particle determine the dynamics of the 
fluid. Assuming that the condensate does not interact with any other form of energy the above pressure 
leads, via the conservation balance, to 

K Po 1.266 X Qbeo X (m^/lmeV)-3 x (/g/lO^fm) 

^''^ 27Th-^la - Po ^° 1 + 1.266 X f}5,o X (mx/lmeV)-3 X (Za/lQOfm)' ^' 

The current value of the scale factor a is ao = 1 and the current fractional density of the BEC dark matter 
is rifceo = Pbeo / Pc wherc pc the critical density. Combining the above relations we obtain the equation of 
state parameter of the BEC dark matter 

Po 

Wbe = ■ (7) 

a-^ - Po 
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A crucial quantity in this model is the moment z^r at which the condensation takes place (i.e. the 
transition from the "normal" DM phase to the Bose-Einstein state). Note that before the transition the 
bosonic DM particles have decoupled from the primordial plasma and have formed an isotropic gas in 
thermal equilibrium. From kinetic theory the pressure of a non-relativistic gas in this regime is given by 

mis] 



a 



2 2 

p C 



/(p)d3p « 47r 



g 



■dp ^ Py 



2 2 



(8) 



where g is the number of spin degrees of freedom, h is the Planck constant, p is the momentum of the 
particle that has energy E — with distribution function /. For the velocity dispersion 

(T^ = (tJ-^) /3c^ we assume a mean velocity square (y^) = 81 x 10^'*cm^/s^, leading to ct^ = 3 x 10~^ [16]. 
By imposing the continuity of the pressure at Zcr as the thermodynamical condition to be satisfied at the 
critical redshift, we obtain [16 

■ (9) 



(1-Po) 



The quantity plays the role of the dark matter equation of state parameter for z > z^r- However, the 
usual approach in cosmology considers "normal" DM as a standard pressureless fluid. 

Since we will assume positive scattering lengths (0 < po < 1) the equation of state parameter Wbe 
can be negative if < po. This can occur in the past. However, typical values of the free parameters 
produces p < 10~^. This means that < for z > 215. On the other hand, p < 10~^ also implies that 
Zcr ~ 2. This means that once the condensate takes place its equation of state is always positive. 

Assuming that the transition occurs during the matter dominated phase the radiation fluid can be 
ignored. The cosmic background expansion H — a/ a after Zcr is given by 



Hi 



^^6eo(l - Po) 



PO 



(10) 



where ffp is the Hubble constant. Assuming the WMAP7 results [33], the current fractional density 
of the baryonic component is f^bo — 0.045 and the density parameter of the cosmological constant is 
57a — 0.73. We will assume that the transition occurs suddenly at the redshift Zcr- This means that the 
effective expansion before such time is governed by the standard cosmology where po = 0. Actually, this 
assumption is not exactly true. The first order transition occurs at a fixed temperature but it takes some 
time At to fully convert normal dark matter into the BEG state. However, as shown in [16j for typical 
values of the BEG parameters the full conversion takes At ^ Wr^Gyrs. This value can be relevant at 
high redshifts but as we will show, concerning the matter power spectrum, the present model displays 
more appreciable features if the transition happens at low redshift where At ^ Vi~^Gyrs is negligible. 



2.2 Perturbations using the neo-Newtonian cosmology 

At scales larger than the horizon the Newtonian theory fails and the full relativistic equations should be 
adopted. However, cosmology can be understood within the Newtonian framework if we properly interpret 
the physical quantities like velocity, energy density and gravitational potential. Such approach is known 
as Newtonian cosmology that is a useful aproximation for a Einstein-de Sitter Universe. However, if the 
inertial effects of the pressure become relevant as for example during the radiation phase or at the onset 
of the accelerated expansion, Newtonian Gosmology also fails and we need a more appropriate set of 
equations. The inclusion of pressure in the Newtonian cosmology in order to make it relevant for the 
homogeneous and isotropic background gave rise to the neo-Newtonian cosmology |24) . 

A relativistic component sources the gravitational field through the energy-momentum tensor 



T^"^ = {p + P)u''u'' + Pg'"'. 
Gontracting T^'^ with and /i^q — gp,a + Ufj,Ua we obtain, respectively 

p + Vr(pw) + PVrW = 0, 



V + (v.Vr)v 



P + P 



Pv 



(11) 



(12) 



(13) 
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Since in the neo-Newtonian cosmology it is assumed that the effective energy of the fluid is the trace of 
T^^ , then the pressure comes into play. Moreover, the gravitational interaction occurs via the modified 
Poisson equation 

VV = 47rG'(/9 + 3P). (14) 

Equations (jl2m4p were used in [17 [ in order to study cosmological perturbations. However, it is ex- 
pedient to expose a cautionary remark. The neo-Newtonian first-order perturbation dynamics and its 
relativistic counterpart coincide in the case of a vanishing sound speed only (251 . Hence, one may expect 
that the correct relativistic results are reproduced by the neo-Newtonian perturbation dynamics on all 
perturbation scales at least for small values of the sound speed. 

The cosmological constant does not fluctuate in the standard approach. Both the BEC dark matter 
and the baryons obey separately the equations (|12p and (fT3|) . The r.h.s of the Poisson equation will 
include the contributions from all components. Our set of equations will be composed of five equations. 
We introduce perturbations in such equations writing each quantity / — {p,v,P,(j)} as / — >■ / -I- f{r,t) 
where the symbol " ~ " (hat) means a first order quantity. This allows to calculate the matter density 
contrast S = p/ p that will be used to compute the power spectrum. 

Collecting the first order terms and Fourier transforming the perturbations as f{r,t) = Sf(t)e~^, 
with k being the wavenumber of the perturbation, we end up with (details of such calculations can be 
found in ^) 



3Wbe 



Ha a? 



\ H a 1 + Wbe 



(1 + 3cf)(l + Wbe) Sbe = 7TT7^(1 + C^)5b, 



where k^^ = 3000/iA/pc is the present Hubble length. In equations (fTSj) and (fT6| the symbol / means 
derivative w.r.t. the scale factor. The density contrast of the baryonic fluid and the BEC dark matter 
are, respectively 6b and 6be- The speed of sound of the condensate fluid is 

2 9pbe „ 

which is a fundamental quantity for the power spectrum. If the BEC behaves as standard cold dark 
matter, equations (jl5p and (jl6p admit the usual solution 6b ^ 6be ^ o,. In the next section the above set 
of equations will be solved for different values of the parameters la and m^. 

3 The matter power spectrum 

The matter (baryonic) power spectrum is defined as 

Pik) = \6b{z = 0;k)\\ (18) 

where 6b{k) is baryonic density contrast calculated in equations (|15p and (|16p at the present time. The 
baryonic agglomeration 6b is driven by the gravitational field which is sourced by all the forms of energy. 

In order to solve equations (|15|) and ()16|) we need to set the initial conditions for 6b and 6be and their 
derivatives at Zcr where the condensate takes place. Since Zcr = Zcrilai'm-x) ^'^^ each chosen couple of 
values {lajfTi^) we need different initial conditions. 

Assuming a primordial Harrison-Zeldovich power spectrum P ^ k and evolving it with the appropriate 
growth function and the BBKS transfer function 27 we obtain the theoretical expression for the power 
spectrum today. We will denote such result by Phz{z)- As the ACDM model also fits the data we 
then integrate back in time the perturbed ACDM equations from z = (where the spectrum Phz{z) is 
assumed) until z = z^r- With this approach we have the standard power spectrum at z^r that will be 
used as initial conditions for the BEC model. Now, we are able to evolve equations and (|16p from 
Zcr until today. 



4 



15 (meV) 
20 (meV) 
25 (meV) 
30 (meV) 
35 (meV) 




Log,ok/h (Mpc)"^ 



■c 
of 



m^ = 300(meV) 

- = 400 (ineV) 

= 500 (meV) 

- ™^ = 600 (meV) 

-my= 700 (meV) 




m^=lS (meV) 
= 20 (meV) 
»i^=2S(meV) 
= 30 (meV) 
m,= 35(meV) - 




0.06 0.08 0.10 0. 



k/h (Mpc)"^ 









































-1.60 -1.55 -1.50 




-1.0 -0.fi 



0.04 0.06 0.08 0.10 0.12 0.14 0.16 



Logiok/h (Mpc) ^ 



k/h (Mpc)-^ 



Figure 1: Matter Power Spectrum. Red line is the Harrison- Zeldovich spectrum with BBKS transfer 
function. BEC cosmology with la = lO^fm {la = 10^°fm) was adopted with different masses in the 
top-left (bottom-left) panel. Top-right (bottom-right) panel shows the relative difference between BEC 
models with la = lO^fm {la = lO^'^fm) and the standard Harrison- Zeldovich power spectrum for various 
masses m^. 



Since the linear scales probed by the 2dFGRS [21] data correspond to k/h < 0.185Mpc ^, the 
final spectrum for the BEC cosmology is normalized in such way that PBEc{kh^^ = O.lS^Mpc"^) = 
PHz{kh-'^ =0.185Mpc-^). 

In the left panels of figure 1 we show the power spectrum for some model parameters. The data 
points come from the 2dFGRS project [21]. In the top-left (bottom- left) panel we have fixed la — lO^fm 
{la = 10^°fm) and calculated the power spectrum for the BEC model PsEcik) assuming different values 
of the mass m^. These values adopted here for la are the usual scattering lengths found in laboratory. 
The red solid line in both panels corresponds to the power spectrum PHz{k) obtained from the BBKS 
transfer function. The right panels show the relative difference {Pbec {k) — Phz {k)) /Phz {k). 

A remarkable point of the present model is that if la = lO^fm {la = 10^°fm) the transition occurs in 
the future, i.e. z < 0, for masses < 5meV (m^ < lOOmeV). This sets a lower bound to m^. that can 
be tested using current observations. 
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The BEC cosmology exhibits an excess of power for any values of the model parameters. The chosen 
masses correspond to the largest deviations from Phz{z) we have found. Although the small differences 
they are indeed perceptible. For example, if la = lO^fm, a typical value for the scattering length, the 
difference can achieve ^ 8% if = 20meV (see top- right panel in figure 1). Deviations of the same 
order have also been found for a mass = 400meV if la = 10^'^fm which is the largest scattering length 
used in laboratory. For the parameters values used in figure 1 we show the corresponding critical rcdshift 
Zcr in table 1. Note that the masses adopted here obey the constraint < 19eV found in Ref. [19) . 







la 


= 10«fm 






la 


- 10^" 


:m 




m^(meV) 


15 


20 


25 


30 


35 


300 


400 


500 


600 


700 


Zcr 


2.17 


3.23 


4.29 


5.35 


6.40 


1.95 


2.93 


3.91 


4.89 


5.88 



Table 1: Values of the masses used in Figure 1 with the corresponding critical redshift Zcr at which 
the transition to the BEC phase occurs. 



4 Conclusions 

Wc have assumed that dark matter is composed of scalar particles that are able to form a Bose-Einstein 
condensate at some critical redshift Zcr- At this point, a first order phase transition drives the conversion 
of "normal" dark matter into the BEC. After Zcr the dynamics of the dark matter component evolves in 
a slightly different way [16^. Consequently, one can also expect a different perturbative dynamics. Indeed 
BEC dark matter accelerates the gravitational clustering at small scales [XT', TS^ but there are claims that 
at large scales the BEC dark matter behaves effectively as a typical cold dark matter component |19j . 

Using the matter power spectrum we have shown that if such phase transition has occurred in the 
recent Universe this process was able to leave small, but perceptible, imprints on the large scale structure. 
Assuming la — lO^fm the BEC dark matter models shows corrections of the order of a few percents for 
masses 15 — 35meV. Adopting la = 10^" fm corrections of the same order are obtained for masses 
300 - 700meV. 

For the relevant parameter values studied here (see table 1) the transition to the BEC phase occurs 
at low redshifts. Since the standard cosmology remains unchanged before Zcr the CMB physics at the 
last scattering surface will be the same. However, the BEC dark matter would modify the gravitational 
potential just after Zcr while the speed of sound is nonzero leading to a contribution to the integrated 
Sachs- Wolfe effect. Such analysis needs the relativistic perturbation theory which is beyond the scope of 
this work. 

Although the small infiuence of the BEC phase on the matter power spectrum a more quantitative 
analysis can be performed using Bayesian statistics. This would estimate the preferred values of the 
model parameters. We leave this analysis for a future work. Also, for a more general study of BEC 
systems we can consider the case e ^ which would provide a richer dynamics for the BEC system. 
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